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ABSTRACT 



Context. Fossil systems are defined to be X-ray bright galaxy groups (or clusters) with a two-magnitude difference between their two 
brightest galaxies within half the projected virial radius, and represent an interesting extreme of the population of galaxy agglomera- 
tions. However, the physical conditions and processes leading to their formation are still poorly constrained. 

Aims. We compare the outskirts of fossil systems with that of normal groups to understand whether environmental conditions play a 
significant role in their formation. We study the groups of galaxies in both, numerical simulations and observations. 
Methods. We use a variety of statistical tools including the spatial cross-correlation function and the local density parameter A5 to 
probe differences in the density and structure of the environments of 'normal' and 'fossil' systems in the Millennium simulation. 
Results. We find that the number density of galaxies surrounding fossil systems evolves from greater than that observed around 
normal systems at z=0.69, to lower than the normal systems by z=0. Both fossil and normal systems exhibit an increment in their oth- 
erwise radially declining local density measure (A5) at distances of order 2.5 r vir from the system centre. We show that this increment 
is more noticeable for fossil systems than normal systems and demonstrate that this difference is linked to the earlier formation epoch 
of fossil groups. Despite the importance of the assembly time, we show that the environment is different for fossil and non-fossil 
systems with similar masses and formation times along their evolution. We also confirm that the physical characteristics identified in 
the Millennium simulation can also be detected in SDSS observations. 

Conclusions. Our results confirm the commonly held belief that fossil systems assembled earlier than normal systems but also show 
that the surroundings of fossil groups could be responsible for the formation of their large magnitude gap. 

Key words. Methods statistical- Galaxies:clusters:general- Galaxies:evolution 



1. Introduction 

The hierarchical structure formation paradigm is successful in 
predicting many of the properties of galaxy groups and clusters. 
However, to-date, the specific conditions and processes that lead 
to the formation of a 'fossil' system have still to be unequivo- 
cally identified. 

'Fossil' systems are defined as spatially extended X-ray 
sources with an X-ray luminosity Lx > 10 42 fc2 erg s _1 whose 
optical counterpart is a bound system of galaxies with AMu > 2 
mag, where AM12 is the difference in absolute magnitude in R- 
band between the brightest and the second brightest galaxies lo- 
cated within half t he projected virial radius of the systems. (r V! >) 
(I Jones et al. 2003). Broadly speaking, this means that fossil sys- 
tems consist of a relatively isolated, luminous, early-type galaxy 
embedded in a swarm of much smaller galaxies and an extended 
X-ray halo. These systems may therefore be of considerable im- 
portance as the place of formation of giant elliptical galaxies. 

The question then naturally arises: how do these systems 
form? One suggested scenario is that fossil systems form when 
a group or cluster remains undisturbed for a significant fraction 
of a Hubble time. Within the context of hierarchical structure 
formation, this means that fossil systems would have assembled 
their dark matter halos earlier than normal groups, thus leaving 



enough time for L* objects to merge into the central galaxy by 
dynamical friction. In other words, any large, centrally located 
satellite galaxies would have been tidally stripped, disrupted, 
and finally cannibalised by the brightest galaxy, naturally pro- 
ducing the characteristic magnitude gap (> 2mag). In the present 
epoch, the central regions of fossil systems would therefore ex- 
hibit the observed lack of L* galaxies, but the fainter end of the 
luminosity function would remain intact, because of the longer 
timescale of dynamical friction for low mass galaxies. 

The next question to arise is then: which systems are 
most likely to experience this kind of evolution? Much ef- 
fort has been devoted to answering this question. Beside 
minor differences in the definition of fossil groups, the 
fossil group phenomenon has been studied from the broad 
range of observational, a nalytical, numerical, and semi- 
analy tica l points of vie w ( Vikhlinin et al ] Il999t Jone s et al 
1 2003 
2006 



D'Onehia etalT. I2005T 



Cvpriano et al 
Milosavlieyic et alJ I2OO6 



Sales et al 



120071: 



iMendes de Oliveira et al. 

2006 



20061: iKhosroshahi et alJ 



van den Bosch et alJ 



Mendes de Oliveira & et al 



2007 



2008 



von Benda-Bec kmann et al.1 

2008h . Most of these works have 



been motivated by the questions of whether the large magnitude 
difference in fossil groups implies that they are a distinct class of 
objects or they simply represent a tail of the cluster distribution. 
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Since the number of fossil systems identified observationally is 
rather low, N-body cosmological simulations are of considerable 
importance, as they allow us to perform statistical analysis. The 
largest cosmological numerical simulation presently available 
is the "Millennium simulation" dSpringel et all 120051 hereafter 
MS). When combined with semi-analytical models of galaxy 
formation, this simulation constitutes a useful tool in addressing 
open issues surrounding the formation of fossil systems. Three 
important studies used this to ol to analyse the e volution of 
fossil groups. On the one hand. iDariush et ail d2007l) concluded 
that fossil systems identified in the MS assembled a larger 
fraction of their masses at higher redshifts than non-fossil 
groups. Therefore, they suggest that the most likely scenario 
for fossil groups is that they are not a distinct class of object 
but simply e x ample s of groups/clusters that collapsed early. 
IDariush et all (1201 Oh then suggested refinements to the fossil 
definition to enha nce its efficiency in detect ing old systems. 
On the other hand, iDfaz-Gimenez et all d2008l) studied the first 
ranked galaxies in the MS, comparing the merger history of 
the central galaxies in fossil and normal systems. By analysing 
central galaxies with equal stellar mass distributions, they found 
that, despite the earlier assembly time of fossil systems, first 
ranked galaxies in fossil groups assembled half of their final 
mass and experienced their last major merger later than their 
non-fossil counterparts, which implies that they followed a 
different evolutionary pathway. Consensus has clearly yet to be 
reached regarding the nature of fossil systems. 

Few works have studied the environment of fossil sys- 
tems or the influence that it could have had on their forma- 
tion, mainly because there are too few known fossil groups 
to allow a reliable statistical analysis. However, observa- 
tionally, it has been suggested that fossil systems inhabit 



Vial chins .-W. v] , 



Matching M v[r and z a 



under-dense regions ( Jones et al. 120031) . On the other hand 



Ivon Benda-Beckmann et all (120081) used numerical simulations 
to conclude that there is no tendency for fossil systems to be 
preferentially located in low density environments. They also 
showed that many galaxy groups and clusters may undergo a fos- 
sil phase in their lives, but may not necessarily still be fossil sys- 
tems at z — because of the infall of L* galaxies from the large- 
scale environment. However, we note that their study was limited 
to groups of relatively low virial masses (1-5 x lO 13 ^ 1 Mo), 
which might not encompass the charact eristic X-ray emissio n of 
observationally selected fossil systems (IDariush et al.ll2007l) . In 
addition, several works on galaxy groups have shown that their 
virial masses and the assembly times strongly depend on their 
environment dGao et alJl2005t iBerlind et aljl2006t IWetzel et all 
l2007tlJing et al.ll2007h . 



In this work, we attempt to obtain a clearer picture of fossil 
systems and their environments by analysing the surroundings of 
systems extracted from the MS. Throughout, we refer to groups 
and clusters that meet the 'fossil criteria above as 'fossil sys- 
tems' or simply 'fossils' and those that do not meet the criteria 
as 'non-fossils' or 'normal systems. 

The layout of this paper is as follows: In Sect. [2] we de- 
scribe the construction of the sample of fossil and non-fossil 
groups. In Sect. [3] we study the outskirts of these groups us- 
ing the two-point cross-correlation function, while in Sect.lU we 
describe their environment by means of the local density profile. 
The analysis of the differences between the environments of fos- 
sils and non-fossils are carried out in Sect. [5] A comparison with 
observationally selected fossil groups is performed in Sect. |6l 
while we summarise our conclusions in Sect. [71 
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Fig. 1. Upper panels show the virial mass distributions of the 
samples of fossil (empty) and AM12 < Q.Smag (filled) groups 
in the MS, while lower panels show the assembly time distri- 
butions. Left panels: "Mass-matched" samples. Non-fossils that 
were selected to reproduce the same virial mass distribution 
of fossil groups in order to avoid a mass bias. Right panels: 
"Assembly-matched" samples. Fossils and non-fossils were se- 
lected to have both similar virial mass and assembly time distri- 
butions. 

2. The samples 

Our samples of fossil an d normal groups were extracted from the 
Millennium Simulation dSpringel et alj|2005l) combined with a 
semi-analy tic model of galaxy f ormat ion. These samples were 
defined in iDfaz-Gimenez et all d2008l) . Here, we briefly sum- 
marise their main characteristics: 

- The MS evolved 10 billion (2160 3 ) dark matter particles 
within a periodic box of 500 ft -1 Mpc in a Lambda cold dark 
matter (ACDM) cosmology with parameters: £l m = 0.25, 
Q A = 0.75, cr 8 = 0.9, and h = 0.73. 

- To ext ract galaxies from the MS, we adopted the 
run of iDe Lucia & Blaizotl d2007l) 's semi-analytic model, 
which provides positions, velocities, absolute magnitudes 
(BVRIK), etc. The final output at z = contained ~ 10 7 
galaxies with absolute magnitudes Mr - Slog h < -17.4 and 
stellar masses higher than 3 x 10 8 h At©. 

- Groups were identified at the z=0 output of the simula- 
tion using a friends-of-friends algorithm in real space with 
a linking length of 0.2 of the mean particle density and 
cross-correlated with the D M halos identified in the MS 
dDfaz-Gimenez et aD(l2008h N >. Their merging history was fol- 
lowed back in time to study the evolutionary behaviour of 
galaxies inside and outside the virial radii of the systems. 
Virial radii and positions of the groups at the different out- 
puts used throughout this work were extracted from the pub- 
licly available information of the MS+semi-analytic modeQ. 



1 The Millennium simulation, performed by the Virgo Consortium, is 
available at http://www.mpa-garching.mpg.de/millennium 



Dfaz-Gimenez et al.: Fossil groups in the Millennium simulation 



3 





Fig. 2. Top panels: Group-galaxy cross-correlation functions of the mass-matched samples {^c-gai)- Upper panels: %c-gai m the MS 
as a function of the normalised group-centric distance. Filled triangles {dotted lines) correspond to fossil groups, while open circles 
(dashed lines) correspond to normal groups with a magnitud e gap AM12 sm aller than Q.Smag. Error bars are computed using an 
analytic formula for the Poisson error in the clustering signal dPeacocklll999t) . Lower panels: Ratios of cross-correlation functions 
(Rg) for fossils and groups with AM12 < Q.Smag. Error bars are computed using the usual formula of error propagation. Panels 
from left to right: Different stages of the group evolution (see upper labels). Bottom panels: Local over-density profile of galaxies 
around groups (A5). Upper panels: A5 in the MS as a function of the normalised group-centric distance. Filled triangles {dotted 
lines) correspond to galaxies around fossils groups, while open circles {dashed lines) correspond to galaxies around normal groups 
with AM12 < 0.5 mag. Error bars represent the 35th and the 65th percentiles in the local density distribution at each distance bin. 
Lower panels: Ratios of local density profiles {Ra 5 ) of fossils and groups with AM12 < Q.Smag. Error bars are computed using the 
usual formula of error propagation. Panels from left to right correspond to the stages of the group evolution (see upper labels). 



Specifically, virial radii were estimated using the virial theo- 
rem, with the virial masses and velocity dispersions obtained 
from the MS database. The virial masses used for these esti- 
mates were the masses within the radius where the halo has 
an over-density corresponding to the value at virialisation in 
the top-hat collapse model for the KCDM cosmology. 
- Fossil systems were selected from the sample of systems 
identified at z = to have virial masses higher than 5 x 



10 13 h~ l Mo and absolute magnitude difference between the 
first and second brightest galaxies (AM12) greater than 2 (in 
the R-band) when considering objects within 0.5 r vlr . This 
lower cut-off in virial masses maximise s the probability tha t 
the selected systems are X-ray emitters (iDariush et al.l2007h . 
The final sample comprises 591 fossils. We note that there 
is a smal l difference between this s ample and the one ob- 
tained bv lDiaz-Gimenez et al.l (12.008b . This is due to the lack 
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of some physical properties of some systems in particular 
redshift snapshots used in our analysis. The virial mass dis- 
tribution of fossil systems can be seen in the upper left panel 
of Fig. \l\(empty histogram). 

- Following the procedure of iDfaz-Gimenez et al. we 
define a sample of non-fossils in order to directly com- 
pare with the fossil sample. The non-fossil sample com- 
prises groups with the same lower limit in virial mass, but 
having AM12 < Q.5mag. This sample of non-fossils com- 
prises 1997 systems. It is well-known that cluster forma- 
tion histories depend strongly on the mass of the systems. 
Consequently, differences in the virial mass distributions 
of the fossil and non-fossil samples could introduce bias 
into the results. Therefore, we extracted a sub-sample of 
591 non-fossil groups with the same virial mass distribution 
(Kolmogorov-Smirnov coefficient > 0.98) as the sample of 
fossils (the "mass-matched" samples; see upper left panel 
of Fig. [TJ. The lower left panel of Fig. [JJ shows the assem- 
bly time distributions (time at which groups have assembled 
half of their final virial mass) for the samples of fossil and 
non-fossil (AM12 < 0.5 mag) groups. 

- To test the effects of assembly time on our statistics, we 
also defined samples of fossils and non-fossils with matched 
assembly times (right panels, Fig [TJ. We refer to these as 
the "assembly-matched" samples (although we note that the 
samples are both assembly-time- and mass-matched). 

- We also studied the remaining sample of normal groups with 
0.5 < AM12 < 2, which comprises 6881 groups. We divided 
this sample of normal systems into 15 sub-samples, each of 
them resembling the fossil sample in terms of the number of 
systems and virial mass distributions. 

We note that we are unable to easily assess the effects of 
using a semi-analytic model other than the iDe Lucia & Blaizoi 
(1200% model used in this work. However, given that all sim- 
ilar models that reproduce the galaxy luminosity function rea- 
sonably well, we wou ld expect any effects to be small. We par- 
ticularly note that the IDe Lucia & Blaizol model is one of the 
best at reproducing the bright end of the luminosity function 
dBertone et al.l l2007h . which is of particular importance when 
measuring the AM12 parameter that defines fossil groups. 

We note that our fossil sample is comprised of groups that 
exhibit AM12 > 2 at z=0, and does not include groups that had 
AM12 > 2 at high er redshift, but became 'normal' ag ain by z=0 
(as suggested by Ivon Benda-Beckmann et al. 1 120081 The pres- 
ence of these groups in our normal group samples can only serve 
to reduce any differences between our fossil and non-fossil sam- 
ples and therefore does not affect our qualitative analysis. 

3. The group-galaxy cross-correlation function 

The spatial group-galaxy cross correlation function, %c g {R), is 
defined such that the probability dP of finding a galaxy in vol- 
ume element dV at distance R from the centre of a group is 

dP = f)[l+£ Gg (R)]dV, 

where rj is the mean number density of galaxies in the whole 
MS. The quantity %cg(R) is therefore equivalent to the radially 
averaged number over-density profile. To increase the statistical 
significance of our results, we averaged the values for all groups 
after scaling to take into account the different group sizes. This is 
achieved by normalising the group-centric distances to the virial 
radius of each group (R = r/r V j r ). We refer to everything that 



lies farther than one virial radius from the group centre as the 
"environment" of the system. 

At each stage in our analysis, a comparison was made be- 
tween normal groups having 0.5mag < AM12 < 2mag and those 
with AM12 < 0.5mag. We found no differences between these 
two samples of normal groups at any stage. Therefore, through- 
out this paper we present the comparison of our fossil sample 
with the AM12 < Q.Smag sample. 

To follow the evolution of the environment around groups, 
we chose three different redshift snapshots in the MS: z= 0.00, 
0.36, and 0.69. The highest redshift value (0.69) represents the 
time by which 90% of fossil groups have assembled half of their 
final virial mass, while the middle redshift value (0.36) repre- 
sents the time by which 90% of non-fossil groups have assem- 
bled half of their final virial mass. These characteristic times are 
inferred from the assembly time distributions shown in the lower 
left panel of Fig. [TJ 

The top panels of Fig. [2] show the cross-correlation func- 
tions obtained for the mass-matched samples of fossil and non- 
fossil (AM12 < Q.Smag) groups, where the evolution of these 
functions is shown by the three redshift snapshots. We explored 
the cross-correlation functions out to a normalised distance of 
~ 10 r V i r which represents a median distance of ~ 17 fr'Mpc 
(see. Sect.|4]i. Since the mean-intergroup separation for systems 
with virial m asses higher than 5 X 10 13 hr x At© is approximately 
20 hr'Mpc dZandivarez et al.l 120031) . this limit in normalised 
clustercentric distances minimises the risk of the spatial super- 
imposition of groups in our main samples. At all stages, all these 
cross-correlation functions show a clear transition from the 'one- 
halo' to the 'two- halo' regimes at the virial radius of the groups 
dYang et alj|2005l) . In general, the one-halo term dominates the 
total contribution in the non-linear regime, while the two-halo 
term captures the large-scale correlations in the linear regime. 

Starting at z — 0.69 (Fig. [31 upper right panel), the number 
density within the groups (< 1 r vi >) is higher for fossils than non- 
fossils. At this redshift, this may be because most fossil groups 
unlike normal groups, have already assembled half of their fi- 
nal masses. Beyond 3 r vir , the environments of fossils exhibit 
higher densities than normal groups. This is consistent with the 
earlier formation time of fossil groups, since halos forming in 
higher density environm ents collapse earlier dSheth & Tormenl 
12004 lHarker et al.ll2006h . There is a transition region between 
1 r V i r and 3 r„> where the environments around fossils and non- 
fossils looks similar (in terms of number density). 

By z=0.36, both fossil and non-fossil systems have clearly 
evolved. The evolution both inside and outside of the virial ra- 
dius is such that fossil and non-fossil systems have become more 
similar. Within the virial radius, this is to be expected since, by 
this time, most normal groups have also assembled more than 
50% of their final mass. 

At z=0, the ratio R^ = ^p g j^ g indicates that fossils appear 
to have slightly fewer galaxies within half a virial radius than 
normal groups. This may simply reflect that the central galaxies 
of fossils have cannibalised their neighbours to produce the two 
magnitude gap. There is again a 'transient' region between 0.5 
r V i r and 2.5 r VJ > where the global densities around both classes 
look very similar. Beyond this region, the environment of fossil 
systems at redshift zero has a slightly lower density than around 
normal systems. 

A possible cause of the local density differences at high red- 
shift is revealed in the lower panel of Fig. [3] which shows the 
evolution of the halo masses with redshift. Although both fossil 
and non-fossil samples reach the same final masses (by construc- 
tion), at earlier times fossil groups can be seen to have higher 
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Fig. 3. Median values of virial radii {upper panel) and masses 
{lower panel) as a function of redshift for fossils {dotted line) 
and AM\2 < Q.Smag groups {dashed line). Error bars are the 
25th and 75th percentiles. 
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4. The local density profile 

To continue our analysis of the environment around galaxy sys- 
tems, we also study the local over-density profile (A(r/r w >)). 
This statistic measures the local density around each galaxy 
and therefore contains information about the local clustering 
(dumpiness) of the groups. 

The local density, rj^, is measured for every galaxy in and 
around groups. It is computed as the number density within a 
sphere of radius defined by the distance {d) to its N-th nearest 
neighbour, i.e. 



3N 
4nd 3 N 



The local over-density is then given by 
&n - — : , 

where fjtf is the mean local density in the simulation. The local 
over-density profile is obtained by computing the median values 
of the over-density associated with each galaxy as a function of 
the 3D-normalised cluster-centric distance, i.e., (A^(r/ r vi >)). The 
N value was tested using two different values of 5 and 10, and 
very similar results were obtained. We therefore show all our re- 
sults using N = 5. As in the previous section, no significant dif- 
ferences were found between the results for normal groups hav- 
ing 0.5 mag < AM12 < 2 mag and those with AM\2 < 0.5 mag. 

The local over-density profile is shown in the bottom pan- 
els of Fig. [2] It can be seen that, as expected, the local over- 
density of galaxies diminishes as we approach and pass through 
the virial radius. However, there is an increase in the local over- 
density around 2-4 r vi > (the "bump"fl This suggests that either 
there is a sudden increase in the global density at this radius, 
or that galaxies at these distances are more clustered than their 
surroundings. The first hypothesis was ruled out by our study of 
the cross-correlation function in the previous section. We note 
that both fossils and non-fossils show this increase in all redshift 
snapshots, although it is more pronounced for regions around 
fossil groups at z = 0. 



4. 1 . The behaviour of the bumps 

In this section, we analyse the position and amplitude of the 
bump observed in the local density profile, focusing on their de- 
pedence on group properties. 



Fig. 4. Local over-density profile of galaxies around fossil {up- 
per panel) and non-fossil {lower panel) groups in the MS as a 
function of the normalised group-centric distance. Filled trian- 
gles correspond to the sub-sample of systems with early assem- 
bly times, while open circles correspond to a sub-sample with 
later assembly times. Error bars are the 35th and the 65th per- 
centiles of the local density distribution at each distance bin. 



masses. This figure therefore suggests that the mass growth rate 
of fossils differs from that of non-fossils. However, these differ- 
ences in halo mass evolution could be a consequence of different 
assembly times, since halos forming earlier are thought to form 
in higher density regions and also to have had more time to gain 
mass. This point will be addressed in Sect. [5] where we discuss 
the mass- and assembly-time matched samples. 



4.1.1. The apparent movement 

When analysing the bumps for the different evolutionary stages, 
it can be seen that they become more pronounced towards z=0, 
but that they also seem to be moving closer to the group. Our 
analysis shows that most of the galaxies located at the position 
of the bump at a given redshift, are also positioned in the region 
of the bumps in later outputs. Thi s finding leads us to rule out the 
possibility of reboun d galaxies dMamon et al.l l2004t iGill et all 
120051 lM amo nl l2006T) being the cause of the bump since these 
galaxies are expected to dramatically change their positions with 
respect to the group centre between redshifts 0.69 and 0.00. The 
apparent movement of the bump could therefore be the result of 
either infalling substructure or simply the growth of the virial 



2 We refer to the region where the local density increases outside the 
virial radius as the "bump". 
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Fig. 5. Same as Fig.|2]but for the sub-samples of fossils and non-fossils with the same virial mass and assembly time distributions 
(samples in right panels of Fig[T]). 



radii of the systems concerned, i.e. the physical location of the 
bump being fixed. The upper panel of Fig. [3] shows the variation 
in the virial radius with time. It can be seen that the median of the 
virial radius for fossils increases from -1.1 Mpch- 1 atz = 0.69 
to -1.7 Mpc /i~'at z - 0.00. The bump in the local density pro- 
file for fossil groups at z=0.69 is positioned at r/r v ; r - 4.5, which 
in physical units corresponds to - 5 h~ 1 Mpc. If the galaxies in 
the bump possess no significant proper motions with respect to 
the cluster centre, then we should find the bump at 2.9 r vi > at red- 
shift z=0.00, i.e., we should see an apparent movement simply 
because the size of the groups increases with time. However, it 
can be seen that the position of the bump at z=0.00 is - 2.5 r vlr , 
which indicates that besides the apparent movement caused by 
the change in size, there is also a real movement of the particles 
in the bump towards the centre of the group. The proper motions 
of galaxies in the bumps have to account for a movement of less 
than 0.5r v ; r , which means that those galaxies have to have radial 



velocities with respect to the centre of the groups of 0.5r w >/Af 
d,50h- l kpc/6M8Gyr = U64kpc/6.U8Gyr - 185 km/s. 



4.1.2. Dependence on the assembly times 

It is well known that both the mass and assembly times of the 
system s are affected by their environment dEspino-Briones et all 
|2007|) . It is widely acc epted that fossil systems assembled 
earlier than non-fossils (D'Onghia et all l2Q05b iDariush et al.l 



l2007t lD faz-Gimenez et al 



2008h . Hence, the following question 



arises: how does the assembly time of the systems affect the lo- 
cal density profiles? To answer this question, we re-sampled the 
mass-matched samples according to their assembly timefl We 



3 According to the definition of iDe Lucia et al 1 l l2006h . we take the 
assembly time z a to be the time when the system has acquired 50% of 
its final (z = 0.0) virial mass. 
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also restricted the samples to Al™ < 10 14 h _1 Ai to avoid differ- 
ences being caused by any differences in mass. Figure [4] shows 
the local density profiles for systems with early assembly times 
(> 80th percentile of the distribution of assembly times)(dotted 
lines), and systems that assembled later (< 20th percentile of z a ) 
(dashed lines). 

The effect of the assembly time on the local density profile is 
clear: Systems that assembled earlier show a more pronounced 
bump than systems that assembled later, and this is true for both 
fossils and non-fossils. 

It can also be seen that the local density profiles of non- 
fossil systems that assembled later (z a < 0.41) show an inner 
peak, very close to the virial radius for the snapshots z=0.36 
and z=0.69 (middle and right lower panels), and this feature 
can also be seen in the local density profile at z=0.69 for fos- 
sil groups with late assembly times (z a < 0.83) (right upper 
panel). The snapshots at z=0.36 and 0.69 represent stages where 
the groups with late assembly time are not yet fully assembled 
(in the sense that they have not reached half of their final mass) 
and are therefore less relaxed. On the other hand, groups (fossils 
or non-fossils) with early assembly times do not show any inner 
peak immediately outside the virial radius. Their earlier assem- 
bly makes these groups more relaxed in the snapshots shown in 
Fig. |4] Therefore, the shape of the over-density profile is then 
influenced by the assembly time of the groups, and also by the 
stage of relaxation of the groups in any given snapshot. 
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Fig. 6. Median values of virial radii (upper panel) and masses 
(lower panel) as a function of redshift for fossils (dotted line) 
and AM[2 < Q.Smag groups (dashed line) with the same mass 
and assembly time distributions. Error bars are the 25th and 75th 
percentiles. 



5. Differences between fossils and non-fossils 

We have demonstrated above that the local density profiles of 
fully assembled galaxy systems increase beyond the virial ra- 
dius. We have found that galaxies remain in this region for a 
long time, gradually approaching the group centres. We have 
also shown that the shape of the local density profile and the 
amplitude of the bumps strongly depend on both the assembly 
times of the groups and their relaxation stage. This last finding 
might lead us to conclude that the environmental differences ob- 
served between fossil and non-fossil systems are purely an evo- 
lutionary effect since fossil groups assembled earlier. However, 
although fossil groups are indeed, on average, older than non- 
fossils, there exist sub-samples of fossils and non-fossils that as- 
sembled at similar times (see lower left panel of Fig. [TJ- Since 
fossils have developed the two-magnitude gap, while non-fossils 
do not, we are led to ask: Could the environment be responsible 
for this magnitude gap in fossils? 

In this section, we therefore investigate whether the envi- 
ronment is different for fossil and non-fossils systems that have 
the same assembly times, and consequently, the same relaxation 
in any given snapshot. To this end, we select new sub-samples 
of fossils and non-fossils having the same z=0 virial mass and 
assembly time distributions (the assembly-matched samples; 
upper right and lower right panels of Fig. [TJ. Both of the new 
samples comprise ~ 330 groups each. 

We first calculate the cross-coiTelation function in order to 
check whether the differences observed in the upper panels 
of Fig. [2] persist. The upper panels of Fig. [5] show the cross- 
correlation function for these new samples. Comparing these fig- 
ures, the strong differences seen, particularly at the earliest snap- 
shot (z=0.69), have become smaller. However, the tendency of 
fossil groups to be in higher number density regions at high z is 
still present, although it is clearly somewhat weaker. To confirm 
this trend, we analysed the evolution of the halo virial masses for 
the assembly-matched samples. In the lower panel of Fig. [6] the 



medians of the halo masses at the different snapshots are shown. 
As can be seen, at z=0, the fossil and non-fossil halos have the 
same mass (by construction). However, at higher z, fossils ex- 
hibit higher halo masses. This is an interesting result since both 
samples have the same final mass and assembled at the same 
time. The result supports the weak trend observed at high z in 
the cross-correlation function. A higher number density around 
fossils clearly allows these groups to reach higher masses than 
non-fossils in the same period of time. Hence, the environment 
of fossils and non fossils are different even for the assembly- 
matched samples, indicating that the differences are a matter not 
only of time evolution. 

We next calculate the local density profiles for the assembly- 
matched samples. The bottom panels of Fig. [5] show the local 
over-density profiles for these new samples. Comparing these 
results with those found in Sect. [4] it can be seen that the lo- 
cal over-density inner peak observed just outside the virial ra- 
dius of non-fossils at z=0.36 and 0.69 is drastically reduced, as 
a result of our having discarded the least relaxed non-fossil sys- 
tems. Moreover, at z=0, the bumps for fossils and non-fossils 
became similar, which indicates that the difference observed for 
the mass-matched samples at this snapshot is due only to differ- 
ences in the assembly times of the groups involved. However, the 
main differences between fossils and non-fossils observed at ear- 
lier evolutionary stages are still present despite there no longer 
being any assembly time dependence. This is a very interesting 
result as it indicates that the environments of fossils and non- 
fossils are indeed intrinsically different at high redshift, and may 
therefore be the cause of the differences in the magnitude gaps 
among these systems. 

6. Comparison with observations 

The most interesting result we have found is related to the lo- 
cal over-density profile at redshift zero. From this profile, we 
observe that the environments of simulated fossil systems look 
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Table 1. Observational Fossil and Normal Groups 



Fossil Groups 


z 


0~r 


• vir 


L x 






[km s~ ] 


[h- 1 Mpc] 


[erg s- [ 10 43 ] 


RXJ 1340+40 18 


0.171 


706 


1.05 


1.24 


p VT1 1 ^Q-l^^ 1 

KAJ 1 1 J7TJJJ 1 


n ns 1 

U.Uo 1 




1 ns 

1 .Uo 


1 1 A 


RXJ 14 16+23 15 


0.138 


500 


1.49 


6.09 


RXJ 1552+20 13 


0.136 


459 


1.22 


2.29 


Normal Groups 










400d-177 


0.169 


680 


1.22 


2.55 


400d-081 


0.064 


449 


1.25 


2.24 


400d-151 


0.160 


502 


1.06 


1.26 


400d-180 


0.167 


452 


1.33 


3.88 



Notes: Normal groups were extracted from the 400d Cluster Catalogue 
http://hea-www.harvard.edu/400d/catalog/table_cat.html 



different from the environments of normal groups. Is it possi- 
ble to observe the same behaviour of the local over-density for a 
sample of observational fossil systems? 

To answer this question, we selected small samples of fos- 
sil and normal systems from a large-scale galaxy survey. We 
used the main spectroscopic sample of galaxies of Sloan Digita l 
Sky Survey Data Release 7 (SPSS DR7,lAbazaiian et all ( 120091) ') 
and, following the work o flVoevodkin et ali i 20091) . selected four 
known fossil groups from their sample. We also constructed a 
sample of four norm al groups from the 400d Cluster Catalogue 
dBurenin e t al. 2007). The normal groups were selected to ensure 
that the groups had similar radial velocity dispersions and mean 
redshifts to those in the fossil sample. The velocity dispersion 
restriction was imposed in order to match the virial mass dis- 
tribution criterion adopted in the simulations, while the redshift 
restriction was designed to reduce the number density depen- 
dence on distance. The virial radii of the systems were computed 
fro m their X-ray luminosi ty following IVikhlinin et al.l (l2009h 
and IVoevodkin et al] d2009h . Details of the systems are given in 
Table [Q 

In the work presented above, we have defined the local over- 
density in the simulation using the full three-dimensional (3D) 
spatial information. This is obviously not possible with obser- 
vational data. In addition, the observable samples of galaxies 
are apparent-magnitude-limited, and the computed local den- 
sities will therefore depend on the redshifts of the groups. 
Consequently, we slightly modified the way of measuring lo- 
cal over-densities in the observational data. Our procedure is as 
follows: 



5. 
6. 



We select a spectroscopic sample of galaxies around each 
group consisting of galaxies within 1000 km/s of the group 
centre to avoid projection effects. 

We then select galaxies brighter than the absolute magnitude 
limit that corresponds to the farthest redshift allowed (v max = 
v r + 1000 km/s) (volume limited sample). 
We measure projected local densities using the third closest 
neighbour (S3). 

We compute the mean projected local density corresponding 
to the redshift of the group £3(2). 

We compute the projected local over-density as A3 



£ 3 fe) 



We compute the projected local over-density profile as a 
function of the normalised projected distance to the centre. 



The projected local over-density profiles in the observational 
data are shown in the left panel of Fig. [7] This figure shows that 
the bump observed beyond the virial radius in the simulated fos- 
sil groups is also present in the observational sample. However, 



it can also be seen that the ratio of the local over-density pro- 
files (Ra 5 ) within 1 r Vi > of systems in the observations differs 
from that shown in Fig. [2] (lower left panel), with galaxies in 
fossil systems exhibiting significantly lower relative local over- 
densities than their counterparts in normal systems. However, 
this can be understood in terms of the differing luminosity lim- 
its applied to the two analyses, i.e., since we consider fossils, 
we expect, by definition, to find a lack of bright galaxies in the 
proximity of the central galaxy compared to similar regions in 
the normal systems, while faint galaxies should be equally dis- 
tributed in both samples. This is only discernible when working 
on observational samples that are magnitude limited, since only 
the brightest galaxies are observed. 

To perform a fair comparison between observations and 
semi-analytic models, it is particularly helpful to construct a 
mock galaxy catalogue. We therefore used the semi-analytical 
model to build a mock sample of galaxies within a cone consist- 
ing of shells constructed from different snapshots correspond- 
ing to the epoch of the look-back time at their distance. Here 
we use the 17 last snapshots, bringing us to a maximum red- 
shift of z — 0.68. For the mock galaxies, we obtain redshifts 
by adding the Hubble flow to the peculiar velocities projected 
along the line-of-sight direction. We compute the observer-frame 
galaxy apparent magnitudes from the rest-frame absolute mag- 
nitudes provided by the semi-analytical model. These appar- 
ent magnitudes are conver ted to the observ er frame using tab- 
ulated k + e corrections (Poggi antl 1 1997b . We set an appar- 
ent magnitude limit R = 17.44, which roughly corresponds to 
the r -sdss band magnitude lim i t of r = 17.77 (see Appendix 
A in iDfaz-Gimenez & Mamonl d2010l) ). We identify groups of 
galaxies in the mock ca talogue following a similar procedure as 
IZandivarez et all d2006l) . We also measure virial masses, virial 
radius, and radial velocity dispersions for the groups following 
the cited work. 

Fossil and non-fossil groups are selected according to the 
criteria described in Sect. [2] We found ten fossil groups in the 
mock catalogue with a median redshift of 0.064 + 0.015 and 
median group radial velocity dispersion of 338 + 78 km s -1 . We 
selected ten non-fossil groups among the 219 identified in the 
mock catalogue to reproduce the group radial velocity disper- 
sion and redshift distributions of the fossil groups. The local 
over-density profile is then measured by using the same pro- 
cedure described for the observational sample of groups and is 
shown in the right panel of Fig. [7] Comparing both panels, it can 
be seen that there is roughly an overall agreement between the 
profiles computed from the observations and the mock galaxy 
catalogue. The only noticeable difference is that the position of 
the bump in the observations is shifted towards larger cluster- 
centric distances. This shifting might be, in part, a consequence 
of analysing slightly different evolutionary stages (z=0.13 vs. 
z=0.06), which would lead to the position of the bump being al- 
tered. However, this difference is also affected by the different 
methods employed to evaluate the virial radii (which are based 
on values estimated from X-ray observations for the SDSS ob- 
servational sample, and the virial theorem radii in the mock cat- 
alogue). 

7. Summary 

We have analysed the environments around fossil and nor- 
mal galaxy systems identified in the Millennium simulation 
dSpringel et al.l 120051) combined with a semi- analytical model 
of galaxy formation ( De Lucia & Blaizotll2007l) . Fossil and nor- 
mal systems have been selected using the same procedure as 
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Fig. 7. Local density profile of galaxies around groups as a function of the normalised group-centric distance. Left upper panel 
show the local density profiles for fossils and normal systems identified observationally (see Table 1) while right upper panel 
show the profiles for Fossils and AMn < O.Smag groups in the MS when introducing a brighter cut-off in the R-band absolute 
magnitude. Error bars are the 35th and the 65th percentiles in the local density distribution at each distance bin. Lower panels show 
the corresponding ratios of local density profiles in each case. 



iDfaz-Gimenez et all (120081) . The samples of normal systems 
were divided into two sub-samples, taking into account the mag- 
nitude difference between the first and second brightest galaxies 
in the systems: AM12 < O.Smag and O.Smag < AM\2 < Imag. 
These sub-samples were constructed to reproduce the virial mass 
distribution of fossils to avoid introducing a mass bias. The re- 
sults obtained for these two samples of non-fossils are almost 
identical. We note that adopting any specific semi-analytical 
model introduces a dependence of the results on the particu- 
lar set of parameters and physical processes that were used in 
the model constructio n. Since this work is intend ed to continue 
the analyses began bv lDfaz-Gimenez et al.l d2008l) . analysing the 
differences caused by using different semi-analytical models are 
be yond the scope of this work . However, we note that the model 
of dDe Lucia & Blaizotl2007l) used in this work was shown to re- 
produce the bright end of the lu minosity function, to which our 
work is most sensitive, very well dBertone et al.ll2007b . 

We have investigated the environments using two different 
approaches, one global, the other local. We firstly studied the 
global density of galaxies as a function of the normalised group- 
centric distance by means of the two-point cross-correlation 
function. Secondly, we focused on the variation in the local 
density around galaxies as a function of the normalised group- 
centric distance by using the local density profile. From the 
global density analysis, we observed that at z=0.36 and earlier 
fossil systems in the MS were surrounded by denser regions 
than normal systems. This result is consistent with the common 
idea that galaxy systems that collapsed earlier lay preferentially 
in higher density regions. Following the evolution of these 
systems to the present day (z = 0), we found that the outskirts 
of galaxy systems in the MS have evolved differentially, finally, 
at z — 0, reaching a point where the regions surrounding fossils 
are slightly under-dense compared to normal systems. This last 
result agrees with some observational resu lts that studied fossil 
groups at low redshifts dAdami et al . 2007b. 



From the local density point of view, the local density pro- 
files have yielded very interesting results. We measured an in- 
crease in the local density profiles in the outskirts of groups for 
both fossil and non-fossil samples at z — 0, but fossils clearly ex- 
hibit a sharper increase. This local density bump is observed for 
fossil systems in every redshift snapshot, decreasing their ampli- 
tude and shifting toward higher distances as we move backwards 
in time. The difference in the local density profiles at z — be- 
tween fossils and non-fossils occurs at ~ 2.5 r/r v ,,>, and indicates 
a clear difference exists between the environmental distributions. 

Even though fossil systems have relatively under-dense re- 
gions in their outskirts according to the global density results, 
the local density analysis implies that galaxy clumps are more 
concentrated around fossils than around normal systems at z=0. 

To address these differences in the environments, we have 
performed several tests: 

- Firstly, in the local density profile, we have found that the 
apparent movement of the bump with time is not due entirely 
to the growth of the groups with time, but also to the motion 
of the galaxies in the over-densities towards the centres of 
the groups. 

- Secondly, we divided the samples of both fossil and normal 
systems into two disjoint sub-samples according to their as- 
sembly times (corresponding to a early and late assembly). 
The local density profiles demonstrate that, for both samples, 
the earlier the groups assembled, the larger the amplitude in 
the local density bump is. Hence, we have demonstrated that 
the formation time of groups is reflected in the amplitude 
of this local density feature. This result is quite interesting 
since the local density profile emerges as a new statistical 
tool suitable for distinguishing the formation epoch of galaxy 
systems at z — 0, and also defining the virialised region of 
systems. 

- Thirdly, to avoid the assembly-time influence on the envi- 
ronment, we have constructed new sub-samples of fossil and 
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non-fossil groups with similar virial mass and assembly time 
distributions. The cross-correlation function and the local 
density profiles for these new sub-samples have shown that, 
even when the groups have formed at similar times, the en- 
vironments around fossils differ intrinsically from those ob- 
served surrounding non-fossil groups at earlier evolutionary 
times. Therefore, this result suggests that environment does 
play an important role in the fossil formation scenario. 

iPorter et"ai] d2008l) studied structures around clusters, and 
found an increment in the star formation activity, peaking at ap- 
proximately twice the virial radius of the clusters. They hypoth- 
esised that an increase in the local density of galaxies at about 
twice the virial radius could explain the increased galaxy star 
formation in terms of galaxy-galaxy harassment. This theory is 
strongly supported by the results found in our work here because 
we have been able to effectively measure the increase in the local 
density that could lead to efficient galaxy interactions. 

Our results also encouraged us to investigate whether these 
features in the local density profile are measurable for observa- 
tional fossil systems. We have computed the projected local den- 
sity profile for four well-known fossil groups and four normal 
groups, both samples having similar physical properties. Our re- 
sults suggest that, even with poor statistics, we can obtain simi- 
lar results to those obtained from the numerical simulations. We 
have found a clear increase in the local density profile for the 
fossil sample. 
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